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Abstract 

The quantum jump approach, where pairs of state vectors follow Stochastic Schroedinger Equa- 
tion (SSE) in order to treat the exact quantum dynamics of two interacting systems, is first de- 
scribed. In this work the non-uniqueness of such stochastic Schroedinger equations is investigated 
to propose strategies to optimize the stochastic paths and reduce statistical fluctuations. In the 
proposed method, called the ^ adaptative noise method\ a specific SSE is obtained for which the 
noise depends explicitly on both the initial state and on the properties of the interaction Hamil- 
tonian. It is also shown that this method can be further improved by introduction of a mean-field 
dynamics. The different optimization procedures are illustrated quantitatively in the case of in- 
teracting spins. A significant reduction of the statistical fiuctuations is obtained. Consequently 
a much smaller number of trajectories is needed to accurately reproduce the exact dynamics as 
compared to the SSE without optimization. 

PACS numbers: 03.65.Yz; 02.70.Ss; 05.10.Gg 
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Recently the description of open quantum systems with Stochastic Schroedinger Equa- 
tion (SSE) has received much attention This is due firstly to the possibility of 
substituting the description of a complex system based on the evolution of its density ma- 
trix by a description based on the stochastic evolution of wave-packets. Therefore these 
methods reduce significantly the number of degrees of freedom to be taken into account. 
A second reason is the increasing computational facilities that allow for massive numerical 
applications. The Monte-Carlo wave-function techniques are presently extensively used to 
treat Markovian master equations in the Lindblad form l|, 0, 0,]^, 0, Q, la S 13 • In that case 
the system density evolution (noted ps) is simulated using SSE on system state- vectors |$). 
Then the dissipative evolution is recovered by averaging symmetric densities p = |$) ($1 
over different stochastic trajectories, i.e. ps = p. Following this strategy, various quantum 
jump approaches based on SSE have been proposed l|, E, Hi- 

Large theoretical efforts have been also devoted to the introduction of non-Markovian 
effects using quantum jumps. A possible way to treat this problem is to use stochastic 
equations that contain a non-local memory kernel ll4J . In that case the average 

evolution of the density matrix is still given by averaging over symmetric densities p = 
1$) ($1, but state vectors evolve according to integro-differential stochastic equations. An 
alternative method that avoids the evaluation of non-local memory kernels is to introduce 
pairs of state- vectors leading to densities of the type p = (^2 1- Then it is possible 
to treat non-Markovian effects approximately using time-local SSE. More recently 
SSE using pairs of state vectors have been proposed as a way to simulate the dynamics of 
interacting systems exactly, with the difference that |\E'i) and |\I^2) are separable states of 
the system and its environment. Such an exact reformulation has already been applied to 
different model cases [la and appears as a candidate to treat the correlated dynamics of 

n n 

self-interacting mesoscopic systems |17l I18| and of systems interacting with an environment 

A common aspect to all stochastic methods used to simulate Markovian and non- 
Markovian dynamics is that the equation of motion of the density matrix does not define 
uniquely the stochastic equation of state vectors (see discussion in ref. ^|). A specific 
choice of equation is generally retained by imposing additional conditions along the stochas- 
tic path. For instance in the Markovian limit, nonlinear state-dependent stochastic equations 
have been obtained by invoking either normalization conditions of state- vectors, orthogonal- 
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jump processes 



0I21L 



or measurement arguments 



22 



2311 ■ In addition the introduction of 



mean-field dynamics jl?! . llS^ significantly improves the stochastic description. Further free- 
dom exists due to the fact that the noise itself can be multiplied by an arbitrary phase-factor 
without breaking the stochastic reformulation. Again, this fact has been noted in different 
studies on Markovian dynamics 111 ] and is sometimes used to pass from one stochastic 



equation to another 



24| . In the non-Markovian case it can also provide appropriate 



modifications of the stochastic dynamics 



25|. 



Although an infinite number of stochastic formulations can generally describe the dynam- 
ics of an open system, the physical interpretation as well as the nature of the process can 
be completely different from one SSE to another. In addition, a series of works dedicated to 
Markovian dynamics with symmetric densities have pointed out that stochastic formulations 
are not equivalent as far as numerical implementation is concerned ^ 18, 26|. For instance 
numerical efficiency can be significantly improved by considering evolution of normalized 
states. The flexibility of SSE approaches has however not been explored yet when pairs of 
states vectors are used. It is clear that, in order to be able to apply the latter theory to a 
large variety of physical problems, specific accurate methods must be developed. 

The main purpose of this paper is to investigate more systematically the freedom in the 
formulation of stochastic process using pairs of states in order to optimize the quantum jumps 
and reduce the number of trajectories necessary to describe the dynamics of interacting 
systems accurately. The paper is organised as follows: first the procedure to reinterpret 
the dynamics of interacting systems in terms of a Monte-Carlo evolution of pairs of wave- 
packets is described. The existence of an infinite number of SSE dynamics leading to a large 
freedom in defining stochastic trajectories is illustrated. In section HTl general strategies using 
this freedom to optimize the quantum jumps and minimize the statistical fluctuations are 
described. In this case we show that the optimal noise depends explicitly on the initial state 
and on the properties of the interaction. In section IIIII we show that the combination of 
mean-field dynamics with the optimization of the noise leads to an additional reduction of 
the statistical fluctuations. Lastly, the method is illustrated quantitatively with a schematic 
model of interacting spin systems. 
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I. SSE FOR THE EXACT DYNAMICS OF INTERACTING SYSTEMS 



Here we give a guideline for the exact reformulation of correlated dynamics in terms of 
diffusive wave-function processes using pairs of wave-packets. Note that a formal derivation 
of the equation given below can be obtained using the Hubbard-Stratonovich transformation 
28 1, see for instance ref. jigj ]. 



A. SSE with pairs of wave-packets 

Following we start from a general Hamiltonian describing the interaction of a system 
and its environment given by: 

Hi = Y.A^®B^ (1) 

a 

Here we assume that the Hamiltonian is possibly already written in the interaction picture. 
Aa and are operators acting respectively on the system and on the environment. Starting 
from an initial uncorrelated state: 

\^) = \<^)(^\X), (2) 

where |$) and \x) are state vectors for the system and the environment respectively, the 
evolution of the system can be replaced by the set of stochastic Schroedinger equations: 

fi|$)=7Ea«a(t)^a|$) 

d\x) =lEaba {t)Ba\x) 

where 7 is a parameter defined formally as 7 = Here a n(t) and ba{t) are complex 

stochastic Gaussian variables that follow the Ito stochastic rules j2^. Note that although we 
consider diffusive Gaussian processes, other stochastic methods could be used like piecewise 
deterministic processes (PDP's) jlfil Isnj]. Indeed, the results presented in this work can be 
adapted to PDP quantum mechanics. In the following, we note X the average of X over 
the stochastic variables. Under the condition: 



aa{tM) = 1, (4) 



the average evolution corresponds to the exact Schroedinger equation 

dt 



d\^lJ) = -Hj\^). (5) 



An interesting aspect of this random process is that stochastic equations Q preserve the 
separabihty of the total state given by eq. ^ so that the procedure can be iterated to 
simulate the evolution on large time scales. 

A second attractive aspect is that the stochastic formulation can be extended to provide 
an exact treatment of the Liouville-von Neumann equation of the density matrix. In that 
case, we consider 



D = \^,){^,\ (6) 

where |\E'i) and |\E'2) is a pair of different system+environment states both given by eq. 0. 
Both states follow independent stochastic equations given (for both) by eq. Q. Noting 



p{t) = D{t), its evolution reads 



^hf^ = [Hi,p] (7) 



which is nothing but the exact Liouville-von Neumann equation. Again, eq. © is preserved 
along each stochastic path. The latter reformulation of exact dynamics is not limited to ini- 
tially uncorrelated states. Indeed if the initial density matrix is correlated, it can be replaced 
by an ensemble of couples of states. The initial density reads: p{t = 0) = (^E'2|) 
where the average here means the average over initial dyadics. Then the complete dynamics 
is obtained by averaging both on the initial ensemble and on the stochastic paths. 

The random process presented in this section describes the exact dynamics of a system 
coupled to an environment. Therefore, although the SSE on wave-functions is Markovian, it 
contains all non-Markovian effects. It is also worth noticing that this method significantly 



differs from quantum jump processes used in the Markovian limit [l 
in the non-Markovian limit with a non-local memory kernel jl^ . ^ 



10| or 



|. The first reason is 

that pairs of state vectors that evolve according to independent SSE should be considered. 
A second important aspect is that quantum jumps applied to Markovian dynamics require 
in general to follow only states of the system. Indeed the effect of the environment has 
already been approximated in the Lindblad equation describing the system density matrix 
evolution. Here wave-packets of both system and environment should be followed in time 
and the system evolution can be obtained using the ps = Tte^p) where Tte^-) denotes the 
partial trace over the environment. However the necessity to follow environment degrees of 
freedom may lead to additional difficulties. 



B. Statistical Fluctuations 



The possibility to perform exact dynamics of interacting systems is of particular interest 
to discuss dissipative effects due to the coupling of a system with its environment. In 



10, hj. However 



particular it does not have the limitations of Lindblad Master equations 
large numerical efforts are needed to treat exactly a physical process due to the number of 
trajectories required to reduce statistical fluctuations of the observables. This is why so far 
it has essentially been applied to rather schematic models. 



Following 



a measure of the increase of statistical fluctuations is given by: 

(8) 



kstat =E[Tr 



D{t)-D it) 



E (Tr {D+ (t) D (t))) - Tr (p^ (t)) 



where \\A\f = Tr {A~^A) is the Hilbert-Schmidt norm. In addition, if the initial state is a 
pure state, Tr (p^ (t)) = 1. 

The growth of statistical fluctuations is of particular importance for numerical implemen- 
tations since it is directly connected to the number of trajectories required to properly repro- 
duce the evolution of a physical system using SSE. Now let us make explicit the evolution of 
the statistical noise associated to the SSE deflned by eq. (jS)). We consider the contribution 



of a single initial couple of state vectors D = (\1/2| and we note Xstat = Tr {D+ (t) D (t)) 
its contribution to Agtat- It is assumed for simplicity that each wave-packet is initially 
normalized. Starting from D, the inflnitesimal increase of Xgtat is given by: 

in which each contribution reads: 

dt 



+2^[a^m{t) (A^)^^ (^^>,J) (10) 



explicitly on the properties of the 
llq speciflc choices of noise have 



with i = 1,2. Therefore statistical fluctuations depenc 
.a„do,„ va.ab.es a. and In seve.a. wo.Us 
been made. In all cases it was possible to demonstrate that the fluctuations have an up- 
per bound that grows exponentially in time. Conjointly, the numerical implementation 
depends strongly on the retained SSE as well as the possible appearance of instabilities in 



calculations 
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26| . Let us consider more systematically the dependence of statistical errors 



with respect to the freedom in the definition of SSE using pairs of state vectors. In the next 
section we specify and use the flexibility in the definition of SSE given by equation Q to 
show that an optimal choice of the noise can reduce statistical fluctuations. 



II. OPTIMAL QUANTUM JUMPS AND THE REDUCTION OF STATISTICAL 
FLUCTUATIONS 

Let us now detail the properties of invariance of the SSE approaches. Starting from a 
general complex noise and that fulfills the necessary condition given by (@)), any new 
couple of random variables {a'^{t), h'^{t)) defined by means of the transformation 

with 

Ca = e''^''^ (12) 

also gives the correct exact dynamics. In the following, Ua and 6^ will be referred to re- 
spectively the scaling factor and the phase factor. Invariance of the stochastic reformulation 
with respect to application of a scaling or a phase shows that an infinite number of SSE's 
exists to simulate the exact dynamics of interacting systems, as it has been shown in most 



of the stochastic theories in Hilbert Space 



m 



l25l |. Such an invariance and in particular 



the invariance with respect to a phase factor has already been noted in several works j^, 12^ 
and has been used for different purposes. 

Let us now apply a scaling and a phase factor on the statistical errors. We see that, while 
statistical fluctuations associated to the initial stochastic equations are given by eq. (jlUp . 
the new SSE leads to different statistical fluctuations given by 



+25? {e^^^'^MWW (Aa)^ {B+)X 



where the indices i = 1,2 have been omitted for simplicity. Thus, the growth of statistical 
fluctuations may significantly differ depending on the noise. In this section we show that 
the scaling factor as well as the phase factor can be properly adjusted to obtain minimal 



statistical fluctuations thus reducing the number of stochastic trajectories in numerical im- 
plementations. It may be seen in the last equation that the parameters of the transformation 
act independently on the two parts of eq. (fT^ and as such, they can be adjusted separately. 
For a component a of the interaction Hamiltonian, the two functions fi^ and Y ^ are defined 
as: 

f^a (Ma) = U^W{t)f UlA^ + —\h^{tt (B^B^ (14) 

r„(0,) = 23? ^e'^'-a^{t)hl{t) (A,)^ (5+) J (15) 

Starting from a given statistical noise [a ait), bait)), the strategy is to find the optimal Ua 
and 6a that minimize respectively these two functions. 



A. Optimal scaling factor 

In eq. (I14|) both {A^A^)^ and {B'^Ba)^ are fixed parameters that depend on the initial 
state. For a given initial state vector, Qa is minimal for 

Under this specific choice the lower limit of Qa reads: 

a 



a = 2\l\aa{t)f ■ \ba{t)f (A+Aa)^ {B+B^)^ (17) 

This expression is valid for any and fulfilling the condition Due to eq. we also 
have \aa{t)\^-\ha{t)\^ > 1. Thus it is better to start with a noise such that \aa{t)\^ ■\ha{t)^ = 1. 
Coming back to the general case, it is important to note that the optimization depends on 
the state on which the expectation values are taken. In particular, if 



the minimal value of Qa corresponds exactly to expression (jl3p . In this case, no reduction 
of the statistical fiuctuation is obtained. Such a case will be presented in a forthcoming 
section. Generally speaking, the more the optimal value of Ua differs from one, the larger 
the reduction of the statistical noise. 
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B. Optimal phase factor 



Let us now turn to the second term in eq. ()13|). In order to minimize the statistical noise, 
the second term should be negative. It is always possible to fix conveniently 9a in order to 
have the minimum value for Fq, i.e: 

(19) 



The absolute value of the expression depends explicitly on 
expression for and ha- We define 

la Xa ~\~ I'Va 



a J) a 



Let us now give an explicit 



(20) 



where all components are real Gaussian stochastic variables. Using condition (0)), we obtain: 

= -1 + 2 (xaX'a + VaV'a + Va^'a^ + XaV'a) (21) 



ajfa 



which leads to the inequality 



(22) 



<-l + 2|a2|.|62|. 

In order to reduce the fiuctuations, expression ()21|) must be maximized. The maximal value 
is obtained if Oq, and 6* are fully correlated, i.e. 6* oc Oq,. Therefore a convenient choice for 
the couple {aa.ha) is 



ha 



(23) 



where = 1 and 5 is a real parameter. The latter expression leads to 



aah*a 



-1 + 2 a? 



\h1\ = 1. Since any noise given by the expression ()23|) can be obtained from an initial real 
noise Xa through the transformation (|TT|l . we simply assume that Then, if 

we define 9ab as: 



the optimal phase factor given by 



20a = 77 



leads to the minimum value 



-2 



(Aa)^ (Bt 



(24) 



(25) 



(26) 
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C. Summary 

In this section starting from the diffusive equation given in section II Al we have shown 
that the noise can be optimized to reduce the statistical fluctuations. For any given initial 
state |vl') = 1$) (g) 1;^), an optimal reduction of the statistical fluctuations corresponds to a 
diffusive process with a specific complex Gaussian noise given by 

- e y "aXo, 
— id 1 




with x\ = \. In these expressions both the scaling factor and the phase factor, given 
respectively by equations (fT^ and depend explicitly on the initial state vector as well 
as on the Hamiltonian. Along the optimal path, the normalization evolution reads: 

dt 



d (xi/ 1 ^) = 2-Y: [^/{AtA^)^ {BtB^)^ - ^{At)^ {A^)^^{Bt)^ {B^)^j (28) 

^4J:F^ (29) 

The procedure described here for a single initial state can be applied to the complete de- 
scription of correlated systems. In this case, an ensemble of wave-packets is considered 
simultaneously, each state evolving according to its own stochastic equation with a noise 
adapted at each time-step. Since for all states the trajectories are optimally chosen, we 
do expect that the total statistical fluctuations obtained by averaging over the ensemble of 
states will also be reduced. In the following, this method will be referred to as the 'adaptative 
noise method'. 



III. NOISE ON TOP OF MEAN-FIELD DYNAMICS 



Before giving an example of application of the adaptative noise procedure, we would like 
to mention that the flexibility on the noise is not the only way to optimize quantum jumps. 
It is possible to introduce a deterministic part in addition to the stochastic contribution 
that partially treats the coupling of the system with the environment without breaking the 
separability of the state vector. One ends up with non-linear, state dependent, stochastic 
equations. Such non-linear equations have been obtained in the Markovian limit by assu miuE 



jumps orthogonal to the state 



20 



21 



or "optimal" measurement arguments 



22, m 
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In all cases we expect that the numerical accuracy will be improved. Closely related is the 



17 



3,Q- 



introduction of mean-field dynamics associated to the interacting system 

In this section the latter method is used. The mean-field associated to the general Hamil- 
tonian given by ^ as well as the associated stochastic dynamics are described. As will be 
seen, the introduction of mean-field gives non-linear equations equivalent to other tech- 
niques. Lastly we discuss how the attractive aspects inherent to mean-field theories and to 
the adaptative noise can be combined to take advantage of both methods. 



A. Optimized Stochastic mean-field dynamics 

Starting from the interaction Hamiltonian the associated mean-field dynamics can be 
obtained using the variational principle 

5{^\ihdt- Hi\^) = (30) 

where is a separable state given by equation Using the variational principle we 
obtain the mean-field equation of motion for each component of the total state vector. 

(31) 

^ \x) = Ea {Aa)^ \x) = hf,p \x) 

where hf.^p and hfjp denotes the mean-field Hamiltonians acting respectively on the system 
and the environment. The mean-field dynamics differs from the exact dynamics since part of 
the coupling is not accounted for. Using mean-field expressions, the complete Hamiltonian 
can be recast as 

Hi I^) = j/^MF ® U + I5 ® hfjp - (A^)^ (B^)^ + E - (^")*) ® - (^-)x)} 1^) 

(32) 

In this expression and 1e correspond to the unity operators acting respectively on the 
system and on the environment spaces. Similarly to the case presented previously, the last 
term in eq. ()32|) can be reinterpreted as an average over stochastic paths, leading to a new 
set of stochastic evolution: 

{§ {hl,p - \ {Aa)^ (5.) J + 7 «a (^a ' (^a) j} |$.) ^^^^ 

d\x)= {§ {hf,p - (Aa)^ (5a) J +lEaba " (5.) J} \x) ■ 
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Thus the introduction of the mean-field prior to a stochastic formulation induces a modifi- 
cation of the operator entering the stochastic contribution. In the latter case Aa and Ba are 
replaced by A'q, = (Aq,— and B'^ = {Ba — {Ba)^)- Similar equations have been obtained 
in ref. using the Hubbard-Stratonovich transformation in conjunction with the Girsanov 
transformation. More generally when mean-field is introduced prior to stochastic formula- 



tion, stochastic terms generally found in the treatment of Markovian 0, 2^ 23|, |30| as 
well as non-Markovian dynamics 3) 3| appear naturally. Such a similitude is not surpris- 
ing since for instance mean-field naturally gives rise to orthogonal jumps. When mean-field 
is introduced, the evolution of the statistical fiuctuations (noted Xgtat), associated to the 
stochastic mean-field dynamics with an initial state |\E'), differs from eq. (fTUI) and reads: 

j\MF _ 



/"Yi \ A. A, / 



It is first observed that the quantum fiuctuations for operators and Ba with respect to 
the state vectors 1$) and \x) appear naturally in this expression. The main advantage of 
mean- field HQ 

is that the latter expression is always much smaller than the first term 
of eq. ()10j) leading generally to smaller statistical fiuctuations. 

Let us now combine the advantage of the mean-field with the optimization proposed in this 
work. At variance with the general case presented previously, we have {A'^)^ = {B'^).^ = 0, 
leading to Ta = 0. As a consequence, only an optimization through the scaling factor can 
be performed. In the stochastic mean-field case, the optimal value is obtained simply by 
replacing respectively Aa and Ba by A'^ and B'^ in expression (fT^. This leads to a reduced 
statistical fiuctuation given by 



dt 



d (VP I ^) = 2-Y: Ji{A+Aa)^ - {A+)^ (Aa)^) [{B+Ba)^ - {B+)^ {Ba)J (35) 
- 2^ E^r- (36) 
Using now the fact that 



{Faf-{F^'^) = {^{A+Aa)^ (5+), (Ba)^ - ^{B+Ba)^ (Aa) ^) > 0, (37) 

we see that an optimized stochastic mean-field dynamics will always further reduce statistical 
fiuctuations. 
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B. Alternative optimization of the phase factor 



In the stochastic mean-field dynamics presented above, there still exists a freedom on the 
phase factor to reduce the statistical errors. More generally although the method proposed 
in chapter ini appears as the natural way to obtain an optimal phase factor, it can obviously 
not be used if = 0. In this section we propose an alternative strategy to obtain the phase 
factor. We assume that cancels out. The optimized statistical fluctuations reduce to 



dt 



(38) 



Here, the SSE without or with mean-fleld are considered indifferently. In the latter case, 
and Ba must be replaced by A'^ and B'^ in the following expressions. For a single time step 
a modiflcation of the phase does not affect directly the expression (jHHjl . However the phase 
factor can be adjusted to act on the value of C2 = {A^Ao,)^, (B^Ba)^ along the trajectory. 
Indeed, starting from the SSE given by equation (jS)), we have ^ 

fEJA;AtA^A^ (B+B^)^ 



(39) 



dc^ 



+ fY.p{AtA^)AB;BtB^Bp 



+i E/3 aphp {AtA^Ap)^ {BtB^Bp) 



■iE.ap^JA^.AtA^lB^B^B^ 



+2f 3? apb*^ {A+A^A^)^ [B^B^B, 

In this expression, only the last term is influenced by a phase-factor. The others are invariant 
with respect to a phase transformation. In order to reduce the statistical fluctuations, it is 
suitable to drive the system towards small values of during the evolution. This can be 



achieved by means of a proper adjustment of 9a- Similarly to section Hi B| we deflne: 



A^A^Ap 



BpB^Bc 



and the new optimal phase factor is then given by 29, 



AtA.AR)^{B}B^B., 



TC 



9'^Q leading to 



^{apb*f,{A:AaAp}^{BlB:B, 



A:^AaA^)^{B^B:Ba 



(40) 



(41) 



As previously we assume that a^b*^ = 1. The method presented here provides an indi- 
rect way to reduce statistical fluctuations when the second term of equation (fTUj) is of no 



^ in the mean-field case, additional terms exist due to the deterministic part in the SSE. These terms are 
not reported here but will not change the conclusion. 
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use. This is always the case in stochastic mean-field dynamics. Therefore in the forthcom- 
ing applications of stochastic mean-field, the phase will always be deduced from the latter 
alternative method. 



IV. ILLUSTRATION 

In this chapter we apply the different optimization procedures described previously on 
an illustrative example consisting of a system of spins in interaction. We show in this ex- 
ample that the statistical fluctuations can be significantly reduced leading to more efficient 
stochastic calculations. A systematic study shows that the reduction of statistical fluctua- 
tions depends on the parameters of the model. 

The analytical model proposed in is considered. The model considers a spin described 
by its Pauli operator a coupled with an environment of spins with spin operators a'^°'^ 
{a = 1, ■ ■ ■ , A^) where is the number of spins. The interaction Hamiltonian is given by 

iJ = 2 5:C„(a+ai")+a_af) (42) 

a 

where Ca is the coupling constant. This model is a simplifled version of the one used for 
the description of a single electron spin in a quantum dot given in jsil. Its simplicity 
is particularly suitable for the present study, enabling us to focus mainly on statistical 
fluctuations. In addition it has a similar form to the one generally taken in the description 



of open quantum systems 



10 



22]. Here, it is assumed that = C/N. In this case starting 



from an initial density p (t = 0) = (0)) (\I' (0)| with 

1^(0)) = !+)® I-,- ■■,-), (43) 

the dynamical evolution is known analytically Q]. In particular the exact system den- 
sity evolution obtained by performing the partial trace on the environment, i.e. ps (t) = 
Tte {p{t)), is given by 

ps{t)=p++\+){+\+p-\-)H (44) 

with p+_|_ = 1 — p = cos^ {2Ct). We consider the model as a benchmark for the different 

SSE considered in this work. In the present case the interaction operators discussed in the 
previous chapters read: Aa = \j (^)cr± and Ba = ^JJ^^^cr^\ 
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We now consider four different classes of simulations. In the first type of calculations 
mean-field is not accounted for. Two calculations will be performed without (noted SSE) 
and with optimization (OSSE). The SSE case can be considered as the reference case and 
corresponds to = = Xa- In the second type of calculations the mean-field is introduced 
and the cases without and with optimization are again considered. They will be referred 
to as Stochastic Mean- field (SMF), again with and Optimized Stochastic 

Mean-Field (OSMF). 

A. Statistical fluctuations 

Let us first consider the case = 1 and C = 0.5. In order to illustrate quantitatively 
the effect of the different optimization, we perform an ensemble of stochastic trajectories 
all starting from the same initial state vector |\E'(0)). The evolution in time of the average 
quantity | \E') is shown in figure [U The SSE, OSSE, SMF and OSMF calculations are 
represented respectively by filled circles, filled triangles, open stars and open squares. In 
each case an ensemble of Ntraj = 10^ trajectories are performed. Let us first focus on the 
stochastic dynamics without optimization. We see that | increases much faster in the 
standard SSE case than in Stochastic mean-field dynamics. As expected the introduction 
of the mean-field alone reduces statistical fluctuations. When optimization procedures are 
introduced, the evolution of (\& | is significantly further reduced. It is worth noticing 
that the OSSE gives a much lower increase than the SMF dynamics. This indicates that the 
effect of optimization on the standard SSE already gives a much better result than the one 
using stochastic mean-field without optimization. Finally Fig. ^ shows that the maximal 
reduction is obtained when both the mean-field and the optimization are taken into account. 

We would like to mention that the present model has several specific aspects. First if the 
phase factor is not optimized (i.e. the phase is set to zero), no improvement is observed. 
Indeed if the original noise is set to Oq, = 6q, = Xa, then (A+Aq,)^ = (B^Ba)^ along the 
trajectory, leading to an optimal scaling factor constantly equal to one. Therefore the strong 
reduction is due to the combined action of the phase factor and the scaling factor. Note 
that in the present model, it is possible to show that for a given initial state 6ab = (^'abi the 
optimal phases are the same in the two optimized calculations. 

In order to quantify the statistical fluctuation growth, each curve is fitted using the func- 
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FIG. 1: Evolution of \ ^) (in log scale) as a function of time t multiplied by the coupling 
constant C for = 1 and C = 0.5 (in arbitrary units). The calculations correspond to SSE (filled 
circles), OSSE (filled triangles), SMF (open stars) and OSMF (open squares). Averages have been 
obtained with 10^ trajectories. 

tion I oc e^°*. A measurement of the increase is directly given by As. We have obtained 
A, - 2.6, 1.3, 0.78 and 0.53 for the SSE, SMF, OSSE and OSMF cases respectively. The 
reduction of the statistical errors has a direct effect on the capacity of the different stochastic 
formulations to efficiently reproduce the exact dynamics: the slower the increase the more 
efficient the stochastic equation. This is illustrated in figure |21 where the evolution of the 
occupation probability n+(t) = {+\ps{'t)\+) is shown as a function of time for the four 
calculations and compared to the exact result. In all cases the same number of trajectories 
Ntraj = lO'' was used. In each case the standard deviation on n+ calculated with the set 
of trajectories is shown in the form of error bars. Large values are observed in the case of 
standard SSE. In addition the SSE calculation is unable to reproduce the exact dynamics 
over the considered time interval. The optimized calculations give a good reproduction of 
the exact dynamics while in the SMF a small departure from the exact solution is observed 
at large time. A precise comparison indicates (as expected) that the OSMF gives the best 
agreement. These results show again that the reduction of the statistical fluctuations im- 
proves the stochastic methods. For a given number of trajectories, the associated standard 
deviations are proportionally reduced. This is illustrated in figure 01 where the standard 
deviations of n+ obtained with the four types of calculations are displayed. Standard de- 
viations follow the same trends as the statistical fluctuations displayed in figure d Note 
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c.t 

FIG. 2: The occupation probability n+ obtained by averaging over stochastic trajectories is shown 
as a function of as a function of time t multiphed by the couphng constant C and compared to 
the exact dynamics (sohd hne) with 10^ trajectories. Top: Results obtained using the SSE with 
(open squares) and without (filled circles) optimization. Bottom, results obtained using the SMF 
techniques with (open squares) and without (filled circles) optimization. In each case the error 
bars represent the standard deviations. 

that all stochastic methods lead to the exact result for an infinite number of trajectories . 
However in order to have the same standard deviations with the standard SSE as the one 
displayed in figure El using OSMF, several millions of trajectories are needed. 

B. Variation of parameters 

We now systematically study the effect of the different optimizations on the behavior 
of statistical fluctuations when the parameters of the model change. In a first series of 
calculations, the effect of the variation of the coupling constant C for a fixed number of 
bath states (A^ = 1) is considered. In the second series of calculations the number of states 
was varied while C = 0.5. was obtained by fitting the statistical fluctuations. We show 
in figure Eland figure the values of As respectively as a function of C and N for the SSE 
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FIG. 3: Evolution of the standard deviations of n+ as a function of time t multiplied by the 
coupling constant C obtained for 10^ trajectories using SSE (filled circles), OSSE (filled triangles), 
SMF (open stars) and OSMF (open squares). 

and the OSMF calculations. In figure lUfor both stochastic calculations, exhibits a linear 
behavior with C in agreement with the expressions (jl(J|) and (j36|) . However the statistical 
fluctuations increase much more rapidly with C in the SSE case than in the OSMF case. 
Therefore it is found that the effect of the optimization increases when the coupling increases. 
Similarly increases faster with the number N of states in the bath in SSE compared to 
the optimized cases. Again the optimization is more efficient when increases. 

In this section we have illustrated by a model case the optimization procedure described 
previously. We have shown that the combined effect of using the mean-field and the adap- 
tative noise significantly improves the quantum jump simulations. Lastly, we have shown 
that the improvement depends explicitly on the coupling strength as well as on the number 
of states in the bath. 



V. CONCLUSION 

In this work the freedom in the stochastic reformulation of the exact dynamics of inter- 
acting systems using pairs of wave-packets has been used to optimize stochastic trajectories. 
For a given initial state an optimal complex noise that minimizes statistical fluctuations is 
obtained. It is shown that the properties of the stochastic variables depend on the specific 
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C (arb. units) 

FIG. 4: Evolution of Xg as a function of the coupling constant C obtained using SSE (filled circles) 
and OSMF (open squares). The result of a fit is displayed with solid (for SSE) and dashed (for 
OSMF) lines and corresponds respectively to Xg oc 5.5C and oc 1.15C 
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Number of states 

FIG. 5: Evolution of as a function of the number of states in the bath for C = 0.5 obtained 
using SSE (filled circles) and OSMF (open squares). The result of a fit is displayed with solid and 
dashed lines for the SSE and OSMF cases and corresponds to A^ oc 1.8\/iV and A^ oc 1.0\/iV. 

state to which it is applied as well as on the interaction properties, hence the denomination 
^ adaptative noise method'' . We have also shown that the method can be combined with the 
introduction of a mean-field prior to the stochastic reformulation. The combination of the 
mean-field and the adaptative noise is expected to lead in general to a further reduction of 
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the statistical fluctuations. 

The different optimization procedures have been studied with a schematic model consist- 
ing of a single spin interacting with a bath of spins. In this model the optimization of jumps 
significantly reduces the statistical fluctuations. Accordingly it has been shown that for a 
given number of trajectories, a much better reproduction of the exact dynamics is obtained 
compared to the standard SSE without optimization. Although promising, more complex 
models should be used to fully demonstrate the power of the method. 

At last we would like to mention that the work described here cannot be directly applied 
to stochastic Manv-Bodv theories to treat self-interacting bosonic or fermionic systems as 
p..oposed .ef. 6 Q. I. th. case, tKe operate. and a. identiea. and a.e 
applied to the same state. Work is currently in progress to modify the original stochastic 
reformulation in order to take advantage of the optimization procedure described here. 

ACKNOWLEDGMENTS 

The author is grateful to Jean-Luc Lecouey, Dominique Durand, Daniel Cussol and John 
Frankland for the careful reading of the manuscript. 



[1] M. B. Plenio and Knight, Rev. Mod. Phys. 70, 101 (1998). 

[2] H.P. Breuer and F. Petruccione, The Theory of Open Quantum Systems (Oxford University 

Press, Oxford, 2002). 
[3] J.T. Stockburger and H. Grabert, Phys. Rev. Lett. 88 (2002) 170407. 
[4] J. Dalibard, Y. Castin and K. Molmer, Phys. rev. Lett. 68 580 (1992). 
[5] R. Dum, P. Zoller and H. Ritsch, Phys. Rev. A45, 4879 (1992). 
[6] N. Gisin and I.C. Percival, J. Math. Phys. A25 5677 (1992). 

[7] H. Carmichael, An Open Systems Approach to Quantum Optics, Lecture Notes in Physics 

(Springer- Verlag, Berlin, 1993). 
[8] Y. Castin and K. Molmer, Phys . Rev. A54 (1996) 5275. 

[9] A. Imamoglu, Phys. Rev. A50, 3650 (1994); A. Imamoglu, Phys. Lett. A191, 425 (1994). 
[10] W. Gardiner and P. Zoller, " Quantum Noise" , Springer- Verlag, Berlin-Heidelberg, 2"*^ Edition 
(2000). 

[11] M. Rigo and N. Gisin, Quantum Semiclass. Opt. 8 (1996) 255. 

20 



[12] L. Diosi and W.T. Strunz, Phys. Lett. A224, 25 (1996). 

[13] L. Diosi, N. Gisin and W.T. Strunz, Phys. Rev. A58, 1699 (1998). 

[14] W. T. Strunz, L. Diosi, and N. Gisin, Phys. Rev. Lett. 82, 1801-1805 (1999). 

[15] H.P. Breuer, B. Kappler, F. Pctruccione, Phys. Rev. A59, 1633 (1999). 

[16] H. P. Breuer, Phys. Rev. A69, 022115 (2004); H.P. Breuer, Eur. Phys. J. D29, 106 (2004). 

[17] O. Juillet and Ph. Chomaz, Phys. Rev. Lett. 88 (2002) 142503. 

[18] L Carusotto, Y. Castin and J. Dahbard, Phys. Rev. A63 (2001) 023606. 

[19] J. Shao, J. Chem. Phys. 120, 5053 (2004). Y. Yan,F. Yang, Y. Liu and J. Shao, Chem. Phys. 

Lett. 395, 216 (2004). 

[20] L. Diosi, Phys. Lett. 112A, 288 (1985). 

[21] L. Diosi, Phys. Lett. 185 A, 5 (1994). 

[22] N. Gisin, Phys. Rev. Lett. 52, 1657-1660 (1984). 

[23] J.K. Breshn, G.J. Milburn and H.M. Wiseman, Phys. Rev. Lett. 74, 4827 (1995). 

[24] G. C. Ghirardi, P. Pearle and A. Rimini, Phys. Rev. A42, 78 (1990). 

[25] A. Bassi, Phys. Rev. A67, 062101 (2003). 

[26] A. Gilchrist, C. W. Gardiner, and P. D. Drummond, Phys. Rev. A 55, 3014-3032 (1997). 

[27] J. Hubbard, Phys. Lett. 3 (1959) 77. 

[28] R.D. Stratonovish, Sov. Phys. Kokl. 2 (1958) 416. 

[29] W. Gardiner, "Handbook of Stochastic Methods", Springer- Verlag, (1985). 

[30] L. Diosi, Phys. Lett. 114A, 451 (1986). 

[31] A.V. Khaetskii, D. Loss and L. Glazman, Phys. Rev. Lett. 88, 186802 (2002). 

[32] C. Cohen-Tannoudji, J. Dupont-Roc, G. Grynberg, Photons and atoms : introduction to 
quantum electrodynamics (New York - Wiley, 1997). 



21 



